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Friedmann-Robertson- Walker universes with a presently large fraction of the 

CN ■ energy density stored in an X-component with wx < —1/3, are considered. We 

J"C , find all the critical points of the system for constant equations of state in that 

i^) range. We consider further several background quantities that can distinguish 

^^ ■ the models with different wx values. Using a simple toy model with a varying 

^sL . equation of state, we show that even a large variation of wx at small redshifts 

i^> is very difficult to observe with di{z) measurements up to z ~ 1. Therefore, 

it will require accurate measurements in the range 1 < z < 2 and independent 
accurate knowledge of O, m ,o (and/or £lx,o) i n order to resolve a variable wx 
i | from a constant wx- 
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1 Introduction 

Since the release of type la Supernovae data independently by two groups, the Super- 
novae Cosmology Project and the High-z Survey Project |T|, 0, indicating that our uni- 
verse is presently accelerating, there has been a wealth of activity trying to explain the 
origin of this acceleration. In particular it was found that, assuming flatness, the best- 
fit universe with a cosmological constant A is given by (fl\,o, ^m,o) — (0-72, 0.28). 
Though these data have still to be confirmed by more accurate ones, as will be the 
case in the near future, it is clear that, if taken at face value, they imply a radi- 
cal departure from usual textbooks Friedman- Robert son- Walker cosmology. Still this 
possibility is very attractive as it seems to make all pieces of our present understanding 
of structure formation fit into a consistent paradigm || ||]. But a true cosmological 
constant, interpreted as a vacuum energy, is about 123 orders of magnitude smaller 
than its "natural" value. This is why there has been a lot of interest in models 
where the dominant fraction of the present energy density is rather some effective 
(slowly varying) cosmological constant term. The prominent candidate in this re- 
spect is some minimally coupled scalar field || || often called quintessence, slowly 
rolling down its potential such that it can have negative pressure. Actually, if this 
candidate is the right one, then it is in principle possible to reconstruct (the relevant 
part of) its potential V(4>) j7|, |3|] and the corresponding equation of state w^ || [TTj 



using luminosity distance measurements in function of redshift z, a procedure that 
can in principle be extended to so-called generalized "quintessence" in the framework 



of scalar-tensor theories of gravity [|11|] . A A-term is equivalent to a perfect isotropic 
fluid with a constant equation of state p = wp, w = — 1. We will be interested in 
this work in some dominant X-component, which can be described as a perfect fluid 
with an equation of state specified by —1 < wx = Px/px < —1/3, thus enabling that 
component to induce accelerated expansion. We compare these models, in particular 
those in the range of values in best agreement with observations, —1 < wx < —0.6, 
with A dominated universes. We don't consider here constant equations of state with 
wx < —1, though this even more mysterious possibility is not excluded by obser- 
vations |T2[ (however we will consider such cases for variable wx)- It is possible to 
have a constant equation of state, even in the framework of quintessence though this 
implies severe constraints on the potential V(<p) f|, [H|. Anyway, results obtained 
here apply to these cases if the X-component is some minimally coupled scalar field. 
Though a constant wx is a simplified assumption, these models give insight into the 
observational properties of more sophisticated models. Finally, we also consider a 
toy model with a variable equation of state where the variation of w(z) is significant 
in the range < z < 2. An arbitrary V{<fi) is expected to give a variable equation 
of state. Our aim is, particularly in view of present experimental and theoretical 
activity, to see whether this equation of state which has non negligible variation in 
the range < z < 2 can be unambiguously separated from models with a constant 
equation of state. We show that the separation of these models using luminosity 
distance measurements will require accurate measurements in the range 1 < z < 2, 
to be probed in the future by the SNAP satellite, a conclusion that can be evidently 
extended to all viable models with variable equation of state. 



2 FRW cosmologies 



Our starting point are the well-known Friedmann equations for a homogeneous and 
isotropic universe, 
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Here the different matter components labelled i, are all isotropic perfect fluids. We 
note from (Q) that a component i can induce accelerated expansion provided pi+3pi < 
0. We are interested first in equations of state of the form p, = Wi Pi, with con- 
stant Wi. It is convenient to write these equations using the dimensionless quantities 
x = — , T = \H \ t and fijn = -^- where H = -, p cr0 = —ft and f stands for 

Oo ' ' ' ' Pcr,0 a,' rl-;,u 87rG j u 

any quantity / evaluated today (at time to), finally fi^o = — 3772-- Note that this 
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parameterization fails when Hq vanishes. Equations (|l|J^) are then rewritten as 

i(^) 2 4* 2 = 5?^^^ + 5"" (3) 
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We see in particular that x = 1, x = 1 or — 1. Further x = — qo, where q is the 
deceleration parameter. Therefore, we have J2i &i,o + ^fc,o = 1 (actually this equality 
can be extended to any given time). We could in principle include the curvature term 
^flk,o into the sum in the r.h.s. of @, it would correspond to Wk = —1/3, but we 
prefer to keep it separate for better physical insight, also this curvature term does 
not appear in (Q). It is suggestive to rewrite (0fj) further as 

-x 2 = E-V{x) (5) 

x = -£, (6) 

ax 

with the following identifications 

v (*) = ~lj2 Q i,o x~ {1+3Wi) (7) 

E = -^fc,o ■ (8) 

In the Hamiltonian formulation of our problem, we have 



where p x is the momentum conjugate to x, with the usual equations of motion 

dH dH 

x = -^— Vx = -- £- • (10) 

op x ox 

Our problem is therefore similar to that of a particle in one dimension with potential 
energy V(x). This gives at once insight into the possible motions, like in problems 
considered in classical mechanics. In the following we consider a universe filled with 
dustlike matter (p = 0) and some unknown component labelled X with negative pres- 
sure, — 1 < wx < —1/3. Hence we have px + 3px < 0. A pure positive cosmological 
constant, i.e. a A-term, corresponds to wx = — 1- 

Present observations indicate flx,o ~ 0.7 and though this is large compared to 
dustlike matter, p x ,o is of course negligible compared to M p ; 4 , the natural order of 
magnitude of the vacuum energy. This is why there are attempts to explain the X- 
component in terms of some scalar field 0, a mechanism quite similar to inflation. 
It is possible in that case too to have scaling solutions with p^ oc x m , m = —3(1 + 
u>0)=constant. However this requires a very particular potential V(<p) for which, 

v(4>) = l — £ . (ii) 

Clearly, when wx < —1/3, the energy density px decreases even slowlier than 
the curvature term oc x~ 2 and, as is the case for a cosmological constant, interesting 
new features appear in the dynamics of our system. We drop in the following the 
relativistic component as it is negligible today and basically of no relevance for the 
issues considered in this work. Hence 



x 2 = n m ,o»" 1 + nx,oa;" (1+3u ' x) + nfc,o (12) 

2 x~ 2" 



x = -~^-X,o(l + 3«,x)x- 2 - 3 ^ , (13) 



and we have readily the identity 



V(x) = ~Q mfi x- 1 - l -n Xfl x~ {1+Zw ^ (14) 



From (13), our universe is presently accelerating provided 



and in particular for a flat universe 

w x < -\n- x ) Q . (16) 

For instance, for Q m fl = 0.3, Qx,o = 0.7, wx < —0.47 is required. Therefore, present 
experimental estimates based on baryons in clusters on one hand giving f2 m0 ~ 
0.3, and on the location of the first acoustic peak in the Cosmological Microwave 



Background detected by the balloon experiments Boomerang and Maxima [O, [D| 
suggesting a nearly flat universe, on the other hand, imply that our universe would 
be presently accelerating for a wide range of constant values, roughly wx < —0.5. 
Accelerated expansion starts at x a given by 

*- 3 l-' = (-l + 3M)§^, (17) 

which corresponds to redshifts 



Za = (-l + 3\w x \)^ -*° -1. (18) 

\"m,0/ 

For — 1 < w — ► —1/3, z a is shifted towards smaller redshifts, for (f2 m> o, ^x,o) = 
(0.3, 0.7), we have z a = 0.671 for a constant A-term, and z a = 0.414 when u> = —0.6. 
The fact that z a is so close to zero, is the cosmic coincidence problem. 

Let us look now for the critical points of our system, points for which x cr = p cr = 0. 
While, by definition, points x a introduced in (|17D have no acceleration, the absence 
of expansion (or contraction) requires in addition 

E - V(x a ) = . (19) 

Any particular point x a satisfying also ([19|) will yield a critical point x cr . Hence, the 
existence of a critical point translates into the following equation for the parameters 

y > o (20) 
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(21) 
(22) 
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a = -3\w x \ (-l + 3\w x \ )^^~ . (23) 

Note that this equation cannot be used in the particular case Q m ,o = 0. For a pure 
A-term, eq. (|20D is cubic. We note first that (3 is bounded from above, — oo < (3 < 1 
(the equality (3 = 1 is unphysical, it corresponds to f2 m)0 = oo). The polynom in ( ]2"DD 
has a minimum in y m i n 

Vmin = ("1 + 3\w X \y^ (24) 



and it is easily verified that, for any given value of (3, or equivalently of f2 m> o, (|20D 
has always at least one solution. We also have the nice equality 

x a = ^ . (25) 

y 

The following cases can be distinguished: 



a) (3 < (0 < Vt m Q < 1): one critical point is found with x cr < 1. 

b) (5 — (O m ,o = 1): there are two critical points. The first one corresponds to 
&>x,o — and the critical point is just x cr = oo. The second one is for y a solution of 
y 3\u>x\ _|_ ay = o an d a; cr < 1. 

c) < j3 < 1 (1 < fi mj o < oo): two critical points can be found with x cr respectively 
larger and smaller than 1. 

d) (3 — 1 (H m ,o = oo): this point is clearly unphysical, however we can consider the 
limit /3 — > 1. In that case, x cr —>■ 1 (for a cosmological constant this solution tends to 
the static Einstein universe). 

An intriguing possibility when there is a constant A-term is a quasi-static stage first 
discussed by Lemaitre |L6| (called "loitering" stage in fL7] ). This stage is followed 



either by recontraction or eternal expansion. Actually such quasi-static stages just 
correspond to configurations in the vicinity of a critical point. Hence such a stage is 
equally possible when — 1 < wx < — 1/3 and it can take place at redshifts 

ztt(-l + 3\w x \)^ y-l = — 1, (26) 

Vmin 

where y is a solution of (|20"D . For Q m ,o < 1, a quasistatic stage in an expanding 
universe is only possible in the past in contrast to expanding universes with Q m ,o > 1, 
where a loitering stage is possible either in the past or in the future, and two such 
universes with f2 m0 = 2 are shown in Fig.[TJ. Let us finally mention that spatially flat 
universes, which are particular subcases of /3 < 0, can have critical points either in the 
infinite past for (fix,o, ^m,o) = (1, 0) or in the infinite future for (fix,o, ^m,o) = (0, 1). 
In the latter case it will take an infinite time for the expansion to relax at x cr = oo 
while in the other case it would take a contracting universe an infinite time to relax at 
x cr = 0. Hence, in a flat expanding universe, and generally in any expanding universe 
with Q m fi < 1, the expansion is eternal for — 1 < wx = constant < —1/3. 

3 Some observational constraints 

It is interesting to see how various models fare when compared with observations. We 
consider in the following several quantities probing the background dynamics that can 
discriminate, and possibly rule out, the different types of universes from each other. 

Age of the Universe: This was actually the main incentive for reconsidering the in- 
troduction of a cosmological constant. Indeed, it is possible for an old universe to 
accomodate with a high value of H . The universe with (f2A,o, ^m,o) = (0.72, 0.28) 
and 0.6 < h = i^o/100 km/s/Mpc < 0.7 has an age 13.7 < to ^ 16 Gyrs, an increase 
of nearly 50% compared to a flat matter-dominated universe! When we change the 
value of Wx, keeping the other cosmological parameters fixed, the age will decrease 
with increasing values of wx- As can be seen from Fig.||], varying for example fix,o 
and wx, we have a degeneracy for the age of the universe for some values of Q m ,o 
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Figure 1: For wx = —2/3 and fi mi o = 2, two different universes having a quasi- 
static (loitering) stage are shown. The two universes have different values Qx,o- The 
present time is chosen here at t — and we show for each universe the relevant 
part, either t < or t > 0, containing a quasi-static stage. The quasi-static stages 
take place in the vicinity of the critical points corresponding to Qx,o — 5.8284... and 
^x,o = 0.17157.... For an expanding universe with il m fi = 2, Qx,o = 5.8284... (plain 
line), loitering is in the past and occurs at z ~ 0.70 while it is in the future at x ~ 3.41 
for an expanding universe with fl m fi = 2, ttx,o = 0.17157... (dashed line). We could 
also have a bounce in the past for slightly higher values of fix,o] the other loitering 
stage can be followed either by expansion (this is the case here, but the expansion 
stage is not shown) or by a recollapse for slightly lower ttx,o values. 




w x 

Figure 2: The age of the universe is shown for various cases in function of wx = —■ 
The central curve corresponds to a flat universe with (fix,o, ^m,o) — (0-72, 0.28). For 
each fix,0y the upper, resp. lower, curve corresponds to r2 m0 = 0.2, resp., 0.4. Note 
that all universes with same fi m ,o have the same age when the X-component scales 
like the curvature term. 



compatible with observations. Therefore the constraint of a minimum age to will 
select a region in the free parameter space (for example the flx,o — % plane if we 
keep Q m fl fixed). We see also the effect of a varying wx when Q m ,o, ^x,o are kept 
constant. Note the particular value wx = — 1/3 for which all universes with same 
flmfl give the same age as can be easily checked with the analytic expression for to- 



te 



dx 



o x H(x) 



(27) 



Of course, we are interested in values of wx sensibly smaller than —1/3. 

Age in function of redshift: this constraint is more detailed than the previous one as it 
puts an age limit at each redshift. This quantity has the following integral expression 



t(z) 



dx 



el x H(x) 



(28) 



An "old" object observed at a certain redshift selects all models with at least that 
age at that given redshift. In that respect, several age constraints have recently 
appeared in the literature JT3]. For example, the age of the radio galaxy 53W091 
observed at a redshift z = 1.55 puts a lower bound of 3.5Gyrs at that redshift. The 
quasar observed at z — 3.62 sets a lower bound of 1.3Gyrs. We see from Fig.[| that 
all universes displayed, besides the flat matter-dominated universe, can accomodate 
these constraints for a value 0.6 < h < 0.7. However a flat purely matter dominated 
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redshift z 

Figure 3: The age of the universe at various redshifts z is shown for several universes. 
The upper curve has Q m ,o = 0.2, the lowest dashed curve (with wx = —0.6) has 
fi m fi = 0.4. The three curves in between all have f2 mj0 = 0.3. All these universes pass 
the constraint set by high-redshift objects mentioned in the text for 0.6 < h < 0.7. 
However the lowest dashed curve, is already on the border for h=0.7. We show for 
comparison t(z) for the Einstein-de Sitter universe (dotted line). 



universe cannot, unless the Hubble constant Hq is lower than the present estimate 
h = 0.65 ± 0.05. If we vary only wx, larger values give lower t(z), hence accurate age 
dating of far objects can possibly set a maximal value for wx- 

The luminosity distance d L (z): it is essentially this quantity (through the measure- 
ment of the Hubble diagram) which was measured for type la Supernovae at redshifts 
up to z ~ 1 by two teams, the Supernovae Cosmology Project and the High-^ Survey 
Project. These two teams have found evidence that our universe is presently accel- 
erating, assuming flatness and wx = w\ = — 1, a best-fit universe is obtained for 
(fl m ,o, ^a,o) = (0.28, 0.72). We have the integral expression 



dr(z) 
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(29) 



with S(u) = sinw for a closed universe, S(u) = sinhw for an open universe and S is 
the identity for a flat universe. In the near future, several groups will gather data for 
samples of around a hundred supernovae each and the SNAP (Supernova Acceleration 
Probe) satellite is planned to make measurements on thousands supernovae with an 
accuracy at the percentage level up to redshifts z — 1.7 |19||. In Fig.HI we see the 



effect of varying the equation of state and/or the balance between the parameters 
^m,0) ^x,o f° r flat universes. Several flat models are considered and some models 
are easily resolved but for small redshifts they are still very close. In particular, 
we see that different combinations of fl m o and wx, can conspire to give practically 
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Figure 4: The luminosity distance (in Hq 1 units) is plotted against redshift for several 
flat universes. For each constant wx value, the upper, resp. lower, curve corresponds 
to fi m fi = 0.2, resp. 0.3. We see that some combinations conspire to give practically 
undistinghishable curves up to redshifts z ~ 1. 



the same di{z) up to redshifts z ~ 1. Therefore high accuracy measurements with 
uncertainties at the percentage level are needed in order to cleanly distinguish the 
models and the need to go to redshifts sensibly higher than 1 is evident. But equally 
crucial, as canalso be seen, one needs independent knowledge of £l m ,o an d ^x,o ( on ly 
one of them if we assume flatness) with uncertainties much smaller than 0.1. To see 
the effect of a variation of the equation of state, we compare several models sharing 
the same value wx{z = 0). As we have said above, if we want to explain the X- 
component using a scalar field, it is natural to expect a time varying wx- Clearly, 
if we are able to pinpoint the equation of state, crucial information will be gained 
about the underlying particle physics model. 

Let us consider a toy model in which a phenomenological wx{z) is assumed and 
we want to consider in how far such a model can be distinguished from a model with 
constant wx- We take the following simple model: 



wx = — 1 + Oi + /3(1 — x) 



(30) 



The two constants a and (3 are chosen arbitrarily. The parameter a measures the 
departure today from an equation of state corresponding to a cosmological constant, 
while the parameter (3 measures the variation in time of wx- Adopting a purely 
phenomenological approach, we will consider also cases where wx becomes smaller 
than —1 which can not be accounted for by a minimally coupled scalar field. The 
energy density px evolves in time according to 



Px = Px,o x 



-3(a+/3) „-3/3(l-x) 



(31) 




redshift z 

Figure 5: The value wx{z) is shown for our toy model. Note that we display here the 
cases with a = which correspond to the present value wx{0) — — 1. 
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Figure 6: The luminosity distance (in Hq 1 units) is shown for some of the variable 
equations of state displayed on the previous Figure and compared to a pure cosmo- 
logical constant {j3 — 0). Flatness is assumed with Q m ,o = 0.3. Obviously accurate 
measurements at redshifts 1 < z < 2, and accurate knowledge of Q m fi are needed in 
order to separate the three cases. 
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From Fig.|G|, it is evident that even a strongly varying vox does not yield a significant 
difference compared to the case with wx=constant. However, at redshifts z ~ (f — 2), 
a measurement with an accuracy at the percentage level might be able to distinguish 
our varying equation of state with (3 = 0.5 from a constant equation of state (and 
same value today). However, for (3 ~ 0.1, clean separation will become a technological 
challenge. 

We see that a variation of the cosmological parameters leads to a better sepa- 
ration of the models than a variation in the equation of state, therefore we need 
independent accurate knowledge of £l m ,o (assuming a flat universe) with an uncer- 
tainty AQ mt Q significantly smaller than 0.1 in order to resolve the various curves. A 
detailed discussion of the feasibility of this resolution, of the use of fitting functions, 
and of the reconstruction of the equation of state, is beyond the scope of this work. 
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